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RECONSTITDTION DE SURFACES EH TROIS DIMENSIOBiS PAR UTILISATION DB 

MQDECfS STATISTIQDES 

La pr^sente invention conceme la reconstitution d' ima- 
ges en trois dimensions et, plus particuliSretnent, la reconstitu- 
tion des cc»itours surfaciques d'une image partir de vues, m&ne 
incoirpldtes, en deux dimensions. 
5 Uki exatple d'e^jplication de la prgsente invention est 

la reconstitutian d' images repr^sentant des 616ments osseux a 
partir d' images en deux dimensions prises, par ex^rple, aux 
raycxis X. De telles images peuvait servir, par ejoaiple, a simuler 
\jne intervention chirurgicale en perraettant au praticien de pr€- 

10 visualiser les positions respect ives des os au niveau d'une arti- 
culation, par exemple, en cas de plastie ligamentaire ou de mise 
en place de prothese- 

Aujourd'hui, pour permettre au praticien de visualiser 
une image en trois dimensions k partir d' images en deux dimen- 

15 sions, on est contraint d'utiliser \m scanner a rayons X effec- 
tuant des tomographies. A partir de ces tonographies, on sait 
reconstituer une image en trois dimensions. Soit il s'agit d'une 
technologie dite de tomographies en trois dimensicais dans 
laquelle un grand nottibre de vues (de I'ordre de 200) sont prises 

20 au moyen d'un scanner rayrais X selon diff§rentes positions, et 
les ccxitours de I'os sont alors reconstituSs pea: traitanent 
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d' linage a partir de ces vues au scanner. Soit il s'agit d'une 
technologie dite de tomographies en deux ditnensions dans laquelle 
un grand nombre de coupes sont prises au moyen d'un scanner a 
rayons X transversalement k I'os, On peut alors reconstituer 
5 l» allure et la structtire de I'os. 

La technique du scanner donne de bons rSsultats ttiais 
est d'lme mise en oeuvre lourde et couteuse. En effet, 1' utilisa- 
tion d'vin scanner pertnet d'dbtenir un ensemble d' images en deux 
dimensions foumissant non seulement des informations sur le 

10 contour mais 6galement sur I'int^rieur de I'os, Or, dans de nom- 
breuses applications, seule la connaissance du contour stirfacique 
de I'os ou de I'cbjet est necessaire. 

Un autre exerrple d ' application de la presente in^v^ention 
est la reconstitution d'ossements incorplets, par exerrple, en 

15 arch§ologie. Une reconstitution d» images en trois dimensions peut 
permettre de retrouver de fa?on quasi-parfaite la forme origi- 
nelle de I'os meme si celui-ci est d^couvert de fagon incoiplete. 
Dans une telle application, des prdbl^s similaires a ceux 
s§s ci-dessus en relation avec la simulation d • interventions chi- 

20 rurgicales se posent. En particulier, il est souvent utile de 
connaitre la forme de I'os sans avoir k se pr§occi^r de sa 

structure interne. 

I»i autre inconv&iient des. techriiques ccainues est 
qu'elles iitposent une dose de radiations itrportante pour le 

25 patient ce qui n'est pas souhaitable. Si cet inconvenient est 
moins sensible dans le domaine de 1 ' arch^ologie oi c'est le cout 
qui prgdomine, il est particuliSrement gSnant dans la simulation 
d' interventions chirurgicales . 

I,a presente invention s» applique plus particulidrenent 

30 a la reconstitution d' images se rapportant a des objets identi- 
fies, c'est-a-dire, dont on connaxt Sl I'avance la forme g&i§rale. 
Par exetiple, pour \m os, il faut au pr§alable decider de quel os 
il s'agit. 

La presente invention vise a proposer un nouveau pro- 
35 c6d6 de reconstitution d' images en trois dimensions qui pallie 
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les inconvgnients des techniques connues. invention vise, en 
particulier, k proposer une solution qui ne necessite pas I'em- 
ploi couteux d'utn scanner a rayons X. 

La pr6sente invention vise egalement a proposer une 
5 solution qui soit conpatible avec \ane exposition tninimale aux 
rayons X ou Equivalents. 

Li 'invention vise en outre a minimiser le nombre de vues 
en deux dimensions nicessaires pour reconstituer 1' image en trois 
dimensions . 

10 Pour atteindre ces dbjets, la prEsente invention pr6- 

voit vn proc6d6 de reconstitution d'une image en trois dimensions 
repr^sentant les contours surfaciques d'au moins un dbjet, k par- 
tir d'au moins une vue en detix dimensions de cet cbjet prise aux 
rayons X, caract6ris6 en ce qu'il consiste k : 

15 determiner la position de la source de prise de vues 

dans un refSrentiel de rfifdrence / 

sglectionxier un module predSf ini constituant une forme 

moyenne de I'objet ; et 

de fagon iterative jusqu't ce que les contours du 
20 module soient tels que les ficarts, entre des rayons de r^tropro- 
jection des contours de l' image en dexxx dimensions depuis la 
source et la surface du moddle, soient minimaux : 

select ionner lone orientation et une position du mqdele 
dans le r^ferentiel de r6f§rence, puis 
25 select ionner ime d6formation du modele pour modifier 

ses contours en trois dimensions. 

Selon un mode de realisation de la pr^sente invention, 
le modele est cbtenu a partir d'une population d'dbjets pour 
laquelle on recherche la correspondance statistique commune & 
30 tous les 6b jets pour determiner vldb forme moyenne et les deforma- 
tions principales par rapport §. cette forme moyenne, de fagon a 
disposer d'au moins un module statistique. 

Selon \m mode de realisation de la presente inventipn, 
les etapes de selection iteratives consistent a faire subir au 
35 modeie statistique, successivement, une transforraaticm rigide 
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modifiant sa position et/ou son orientation et une deformation 
non-rigide modifiant ses contours surfaciques. 

Selon un mode de realisation de la presente invention, 
les contours de 1' image en deixK dimensions sont dbtenus automati- 
5 quCTient en pro jet ant le module dans le plan de 1' image en deux 
dimensions, et en deformant les contours projetSs de fagon a les 
faire coincider avec les points de fort gradient en niveaux de 
gris de 1» image en deux dimensions. 

Selon xm mode de realisation de la presente invention, 
10 la determination automat ique des contours de 1' image en deux 
dimensiCTis est effectuee de fagon iterative, chaque iteration 
etant intercaiee entre deux iterations successives des etapes de 
selection. 

Selon un mode de realisation de la presente invention, 
15 on determine, dans le referent iel de reference, des coordonnees 
en trois dimensions de points de I'objet, de fagon a disposer de 
points de reference suppiementaires pour les etapes de selection 
iteratives de position, orientation, et deformation. 

Selon un mode de realisation de la presente invention, 
20 on utilise plusieurs images en deux dimensions pour lesquelles 
les positions respectives de la source de prise de vues sent tou- 
tes determinees dans le referentiel de reference, et on effectue 

les etc^s de selections iteratives en. tenant cotrpte des rayons 

de retrqprojection des craitours de toutes les images en deux 

25 dimensions. 

Selon un mode de realisation de la presente invention, 
le nombre d' images utilise est fonction de la precisicxi souhai- 

tee. 

Seloa un mode de realisation de la presente invention, 
3 0 la surface du modSle est constituee d' elements de triangle, les- 
dits ecarts etant mesxires par report k des points de certaines 
aretes constituent des generateurs du contour en trois dimen- 
sions. 

Selon iHi mode de realisation de la presente invention, 
35 le precede est applique a la reconstitution des contours surfaci- 
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ques de plusieurs objets li€s entre eux par des relations de 
transformation rigide et/ou Slastique. 

Selon \m mode de realisation de la pr6sente invention, 
le proced6 est appliqu6 a la reconstitution d' images d'os. 
5 L' invention conceme egalement xan systeroe de traitement 

d' images, coctportant des itoyens pour la mise en oeuvre du proc6d6 
de reconstitution d' images en trois dimensions. 

Ces objets, caracteristiques et avantages, ainsi que 
d'autres de la presente invention seront exposes en detail dans 
10 la description suivante de modes de realisation particuliers 
faite a titre non-limitatif en relation avec les figures jointes 
parmi lesquelles : 

la figure 1 repr^sente un mode de realisation d'lm sys- 
teme de reconstitution d' images en trois dimensions k partir d'un 
15 dispositif de prise de vues en deux dimensions selon la pr6sente 
invention ; 

la figure 2A illustre la prise d'une vue en deux dimen- 
sions pour la mise en oeuvre du procSde de reccnstituticn d' ima- 
ges en trois dimensions selon la prSsente invention ; 
20 la figure 2B illustre la prise de trois vues en deux 

dimensions pour la mise en oeuvre du procfide de reconstitution 
selon la prgsente invention ; 

les figures 3A, 3B et 3C illustrent trois etapes carac- 
teristiques du procede selon un mode de realisatiai de la pre- 
25 sente invention ; 

la figure 4 illustre \m mode de determination de gene- 
rateurs de contour selon la presente invention / et 

la figure 5 illustre un mode de determination prefere 
selon 1' invention de I'ecart entre un modele en trois dimensions 
30 et une vue en deux dimensions, 

Les memes elements ont ete designes par les mimes refe- 
rences aux differentes figures qui ne sont pas §. I'echelle. Pour 
des raisons de clarte, seuls les elements qui sent necessaires a 
la conprehension de 1' invention rait ete representes aux figures 
35 et seront decrits par la suite. En particulier, les mpyens de 
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traitement d' images corrprenant les calculateurs informatiques 
n'ont pas et€ representes et ne seront pas d6crits de fa<;on 
detainee, car leur realisation est a la port^e de I'hcmme du 
metier a partir des indications fonctionnelles donnees ci-aprSs, 
5 IMe caractSristique de la prSsente invention est de 

rechercher la position et 1 ' orientation d'un cbjet a reconsti- 
tuer, et dont on connait au moins une vue en deux dimensions, S 
partir d'lme base de donnSes contenant des modules de cet objet. 
Quand plusieurs vues en deux dimensions sont utilisSes, ces vues 

10 sont toutes r^firencees dans vn m&ne r6f6rentiel. Ainsi, 1' inven- 
tion conceme la reconstitution d'une image reprgsentant un objet 
d6ja identifiS et dont on peut disposer de modules de taille 
et/ou de forme diff6rentes. 

une caract^ristique d'xan mode de realisation prefSre de 

15 la presents in:\rention est d'utiliser au moins \xn modele statisti- 
que dSforraable, 6tabli IL partir de la base de donnSes, pour 
reconstituer la forme en trois dimensions de 1' objet. Ainsi, 
1« invention prgvoit de d6finir, avant toute reconstitution, une 
base de donnees contenant des modeles en trois dimensions de 

20 I'dbjet a reconstituer ou, de preference, \m ou plusieurs moddles 
statistiques d^fontables a partir de cette base de donnees. 

La figure 1 represente lane vue schSnatique d'un systSme 
de reconstitution d 'images en trois .diinensi<3ns seslcjn \in mode de 
realisation de la pr^sente invention. L' example de la figure 1 

25 conceme la reconstitution de 1' image d'un os 1 a partir de 
radiographies en deux dimensions. En figure 1, on a repr§sent6 
schgmatiquement et partiellement le corps d'un patient p dont la 
jambe contient I'os 1 que I'on souhaite visualiser. La jambe du 
patient p (done I'os 1) est, par exCTple, pos6 sur une table 2. 

3 0 Utie potence 3 porte un dispositif 4 de localisation en trois 
dimensions qui peut §tre optique, magnitique, mScanique ou ultra- 
sCTiore et qui repere la positicai et 1 ' orientation de multiples 
capteurs-^netteurs CQnstitu6s, par exernple, de diodes infra- 
rouges, de reflecteurs, d'emetteurs magngtiques, a ultrasons, 

35 etc. Ito tel localisateur est parfaitement connu et ne sera pas 
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detaille plus avant- On notera siitplement que, pour garantir une 
localisation correcte du patient (de l*os 1), celui-ci est gene- 
ralCTient ggalement gquipe d'un capteur-emetteur 18 detectable par 
le localisateur 4. En effet, selon 1' invention, toutes les vues 
5 doivent pouvoir etre exploitees dans un meme referentiel (de 
r^fgrence) , qui est associe a I'objet que I'on reconstitue. 

Plusieurs capteurs-^tnetteurs peuvent etre fixes sur le 
systeme de radiologie, proche de la source 7 (capteur-emetteur 
19) ou proche d'un detect eur d'iinage 9 (capteur-emetteur 5), de 

10 fagon a reperer la position du systeme de radiologie par rapport 
au referentiel de reference du capteur-emetteur 18. 

Dans certains cas, les capteurs-&rietteurs 5 et 19 sont 
diff icilement rep^rables en raison de leur eioignement ou de la 
presence d'dbjets parasites dans le chanp de mesure du localisa- 

15 teur 4. Dans ce cas, \m capteur-6netteur 20 peut itre installs 
sur la table 2 dans le chanp de mesure du localisateur 4. On 
positionne suLors le systeme radiologique dans de bonnes condi- 
tions de mesure et cxx repere les positions des capteur-&netteurs 
5 et 19 par rapport au capteur-^netteur 20 une fois pour toutes 

2 0 (on ne renouvelle cette itape que si on bouge le syst&ne de 

radiologie dans son ensemble) . 

Le syst&ne radiologique §tant SquipS de codeurs angu- 
laires sur ses . axes, conwe cela. sera detailie plus loin, les 
changements de position relative du systeme radiologique sont 
25 raesur6s a I'aide de ces codeurs et peuvent ainsi etre reportSs 
dans le r6ferentiel du capteur-emetteixr 20. Pour chaque prise de 
vue, on mesure la relation g§ometrique entre les capteurs- 
emetteurs 20 et 18 et, par ce biais, on reporte la position du 
systeme radiologique dans le rgferentiel de reference du capteur- 

3 0 emetteur 18. Ainsi, toutes les radiograpbies sont calibr§es dans 

le r6f6rentiel du capteur-Sraetteur 18 qui peut 8tre mobile. 

Le dispositif 6 est, par exetrple, constitu6 d'lone 
source 7 a rayons X portSe par une prenidre extr&nite d'un bras 8 
en demi-cercle dont 1« autre e3ctr6mit§ est destinfie k recevoir le 
35 film 9 d'ittpression de la radiographie, ou un capteur elect roni- 



que equivalent tel qu»un anplif icateur de brillance, ou un detec- 
teur plat au silicium ainorphe. La relation entre le bras 8 et la 
table 2 est telle que cette demiere se trouve entre la source 7 
et le capteur 9. Le bras 8 est nonte a rotation autour d'un axe 
5 10, inotorise ou t dgplacettient inanuel, et porte par una potence 
11. L'ensernble de prise de vues en deux dimensions peut ainsi 
toumer autour de I'os 1 pour effectuer le noctibre d€s±r€ de 
radiographies de celui-ci, Le cas 6cdi6ant, le localisateur qpti- 
que 4 peut St re associ6 k un dispositif de codage angulaire de la 

10 position de I'axe 10. 

L'ensennble est pilots par un systdroe infoxmatique, par 
exetiple, ion ordinateur 12 associS, de pr6f6rence, a \m §cran de 
visualisation 13. A la figure 1, on a syttibolisg par des liaisons 
unifilsdLres 14, 15, 16 et 17, les bus informatiques d'Schange de 

15 signaux Slectriques de commande et de donn§es entre 1' ordinateur 
12 et, respectivement, le localisateur 4, la source 7, le tnoteur 
ou le codeur pptidnnel de I'axe 10 et le capteur 9. 

Les figures 2A et 2B illustrent deux nodes de prises de 
vues selon dexax modes de realisation de la presente invention. 

20 Aux figures 2A et 2B, les vues en deux dimensions ont ete repre- 
sentees par les plans respectifs dans lesquelles elles sont pri- 
ses, c'est-a-dire par 1' allure de la surface du capteur 9 (figure 
1) lors des prises, Cela _ correspond aux images radiographiques en ., 
deux dimensions recupgrges par le systeme inf ormatique . 

25 La figure 2A reprSsente le cas d'une seule prise de vue 

dans un plan PI, foumissant une image II de I'os 1. La position 
de la source 7 a §t§ symbolisSe par un point oil convergent des 
rayons r de rStrqprojection des quatre coins du plan PI. 

La figure 2B repr^sente le cas d'xine triple prise de 

30 vues dans des plans PI, P2 et P3, foumissant trois images II, 12 
et 13 en dexix dimensions de l"os 1. En figure 2B, la position de 
la source n'est pas constante, celle-ci est diff^rente pour cha- 
que prise de vue. Toutefois, grice au localisateur 4, toutes les 
positions de la source sont connues dans le referentiel de r6fe- 
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rence. Par soucis de clart§, seuls les rayons rl des quatre coins 
du plan PI dans la position de 1' image II ont §t6 repr§sent6s- 

On notera que la source de prises de vues peut svibir, 
entre deux images en deux dimensions, d'autres tnouvements que 
5 dans un m&ne plan cotrme illustrS par les figures 1 et 2B. En 
d'autres termes, le bras 8 du systeme de prises de vues 6 peut 
presenter plus de deux degr6s de liberty, dont chacun pourra §tre 
TOuni d'un dispositif de codage angulaire. Selon un mode de reali- 
sation pref6r6, le bras 8 est raonte sur deux axes de rotations 

10 horizontal et vertical et un axe de translation vertical. 

On notera ^alement que les images en dexax dimensions 
peuvent n'etre que partielles. Par exemple, dans 1' application 
aux radiographies, celles-ci peuvent etre interpr6t6es par l'qp6- 
rateur pour valider les contours des vues en dexix dimensions d. 

15 prendre en cotipte pour la reconstitution. Une telle interpreta- 
tion n'est pas genante en raison du faible nonibre de vues en deux 
dimensions nScessaire selon 1' invention (moins de dix en gene- 
ral) . La validation des contours dans le syst&ne informatique 
peut §tre effectuSe, par exetnple, au moyen d'lme souris, d'un 

20 crayon optique, d'un Scran tactile, ou equivalent, de fagon clas- 
sique pour \m enregistrement de contour sur une image en niveaxax 
de gris. 

Selon un -autre mode de r&G-isation prefSrS, la determi- 
nation des contours sur les images radiographiques est autcmati- 

25 see en mettant en oeuvre vn precede dit de recalage en deux 
dimensions. Un tel procede consiste a determiner le contour de 
fagon automatique par analyse des gradients les plus forts en 
niveaux de gris. Ce precede est decrit, par exerrple, dans la 
these de Gelu lonescu, presentee publiquement le 4 decembre 1998 

30 a I'lMiversite Joseph Fourier de Grenoble (France), et ayant pour 
titre "Segmentation et recalage d' images echographiques par uti- 
lisation de connaissances physiologiques et morphologiques" . 
Selon 1» invention, on met en oeuvre ce procSde de recalage en 
deixjc dimensions en corabinaison avec vine projection, sur 1' image 

35 en deux dimensions, du modSle statist ique deformable de I'dbjet 
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dont on souhaite reconstituer les contours surfaciques en trois 
dimensions. On effectue une analyse iterative d'lm contour pro- 
jete du module sur 1' image en deux dimensions jusgu'^ dbtenir, 
apres plusieurs passes (par exemple de 3 a 10) , uns identite 
5 entre le contour du modele projete qui sera dgcrit ci-aprSs et 
dont 1' orientation et la deformation auront 6t6 adapt^s par le 
precede de 1' invention, et le ccaitour determine par analyse des 
gradients les plus forts en niveaux de gris. Cette ccnibinaison 
sera mieux cotrprise par la suite en relation avec la figure 5. 

10 Pour siniplifier la pr^sente description, on considirera 

ci-apr^s que les contours des vues en deux dimensiois SCTit 
connus. En pratique, les Stapes du proc6de de 1' invention seront 
effectu^es en boucle, y conpris le receaage en deux dimensions, 
jusqu'S. dbtenir un rSsultat satisf aisant . 

15 Selon 1' invention, une fois que la ou les vues en deux 

dimensions ont et6 dbtenues, il reste a determiner, par recherche 
de correspcndance dans la base de donnees, la forme et la taille 
de I'objet (ici, I'os) en trois dimensions. 

Pour ce faire, 1' invention prevoit de rechercher le 

20 module pour lesquelles les distances sgparant les contours de 
chaque vue en d&jx dimensions et la surface du module scait mini- 
males. 

Ejans le. mode de realisation .prefere qCL on utilise ton 
module statistique deformable (definissant ainsi une famille de 

25 modeies) de I'cbjet k reconstitiier, on. recherche a s'^procher le 
plus possible, par deformation iterative de ce module, de la 
forme 6aat les ccMitours sont tels cjue, cocnme on le verra par la 
suite en relation avec la figure 5, les rayons de retroprojection 
s'c«>puyant sur les points des contours de 1' image (ou des images) 

30 en deux dimensions depuis la source (ou depuis leurs sources 
respectives) sont tous tangentiels a la surface du module, etant 
entendu que les differents contours en deux dimensions sont 
reperables dans le referentiel de reference. 

La realisation du module statistique deformable en lui 

35 meme ne fait pas I'objet de la presente invention et est parfai- 
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tement classique. On se bomera S rappeler que la rialisaticai 
d'un tel module fait g6n6ralement a^Jel k une recherche d'une 
forme moyenne d'une populaticai d'bbjets de minie type (par exem- 
ple, des f^nurs) ccaistituant la beise de donn§es, suivie d'une 
5 analyse en composantes principales pour determiner les deforma- 
tions principales (essentielles) a appliquer a cette forme 
moyenne et obtenir adjisi le module statistique. 

Une mSthode de determination d'un module statistique t 
partir d'une famille d • echantillons est decrite, par exeitple, 

10 dans 1" article "Building a Conplete Surface Model from Sparse 
Data Using Statistical Shape Models: implication to Computer 
Assisted Knee Surgery" de Markus Fleute et Stephane Lavall^e, 
paru dans medical IMAGE (SOMPUTING AND COMPUTER-ASSISTED 
INTERVENTION - MICCAI'98, Springer- Verlag LNCS Series, pages 880- 

15 887, en Octobre 1998, et dont le ccxitenu est incorpori ici par 
reference. On notera que la precision de la reconstruction des 
contours surfaciques operee par !• invention depend des 6chantil- 
lons utilises pour construire la base de donnees du module sta- 
tistique. Par exemple, si un module statistique de femur est cree 

20 a partir d'une populaticai d'echantillons ayant des formes norma- 
les (sans pathologie) , on pourra reconstituer avec precision des 
formes normales mais la precision sera limitee si le femur radio- 
graphie a une forme pathologique . Par centre, si le module statir 
stique est cree a partir d'une grande papulation d'echantillons 

25 contenant a la fois des formes normales et des formes pathologi- 
ques, on pourra reconstituer, avec precision, des dbjets normaux 

et pathologiques . 

Selon 1' invention, une fois que I'on dispose des images 
en deux dimensions, on recherche le cas echeant a quel modeie 

30 contenu dans la base de donnees ces images se rapprochent le plus 
(par exemple, si plusieurs types d'os ou si plusiexars modules 
statistiques d'un mime os sont disponibles) , puis on applique des 
transformations a ce module pour s'approcher le plus possible 
d'vine image en trois dimensions, telle que les rayons de retro- 

35 projection des images jusqu'a leurs sources respectives sent tan- 
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gentiels axxK surfaces de 1' image en trois dimensions, placee sur 
le trajet de ces rayons et dont ils sont respectivement les plus 
proches . 

De preference, on commence par faire subir au moddle 
5 initial au repos, c'est-a-dire correspondant: §l la forme moyenne, 
une transformation dite rigide, c'est-S-dire, ne touchant qu'^ 
son orientatiCTi et sa translation dans l"espace. Puis, quand on a 
dbtenu la posit icxi pour laquelle les rayons de rStroprojection 
sont tous a une distance mdLnimale des contours du module initial, 

10 on fait sxabir \ane deformation ncxi-rigide au modSle, c'est-^-dire 
que sans modifier son orientation, on modif ie sa forme ^ partir 
des donnSes contenues dans la base statistique en modif iant les 
coefficients des modes principaux du module, jusqu'a obtenir la 
forme dont les contours se rapprochent le plus des rayons de 

15 r6troprojection. Le Ccls echeant, le modele est pr§positionn6 de 
fagon trds approximative par I'pp^rateur par rapport aux vues en 
deux dimensions affichees a I'ecran. On notera que les moyens 
(informatiques) de traitement d' images utilises sont classiques 
dans leur structure et n'ont done pas besoin d'etre detailles. 

20 Les figures 3A a 3C illustrent ces deux 6tapes du pro- 

c6de de 1' invention. La figure 3A reprfisente ume vue en trois 
dimensions d'un moddle 21 avant toute deformation, II s'agit, par 
exerrple, d'lm module statistique .^e .1 linv^emtion. positiorme par le 
praticien dans une orientation ^proximative . La figure 3B repre- 

25 sente le module 21' de la figure 3A Sl 1' issue de I'gtape d'orien- 
tatiCTi sans modification de forme. La figure 3C repr^sente 
1» image 22 en trois dimensions resultant de la raise en oeuvre de 
1" invention, c'est-^-dire correspondant au modSle 21', dSforme 
pour que les rayons de r6troprojection soient (idfealement) tan- 

30 gents Sl tous ses contours. 

Par souci de simplification, un seul groupe de rayons r 
de rStrcprojection depuis une source 7 ont et§ represent6s aux 
figures 3A §i 3C et I'ecran de projection de 1' image correspon- 
dante n'a pas 4te represent^. Comme il ressort des figures 3A a 

35 3C, xan grand nonibre de rayons r de la figure 3A ne sont pas tan- 
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gents aux contours du module 21. Le nombre de rayons de r6trqpro- 
jection utilise depend de la precision souhaitee et des caract6- 
ristiques du module statistique, essentiellement, du ncttibre de 
points surfaciques choisis pour servir de reference lors de la 
5 definition de ce noddle statistique. 

De la figure 3 A a la figure 3B, la transformation est 
"rigide", c'est-a-dire que le tnodele n'est pas deforme. Lors de 
cette transfontiation, le moddle 21 subit des translations et 
rotations afin d'obtenir des ttiesiires d'ecarts tnininium entre 

10 chaque rayon de r^tropro jection partant des points de contour de 
1' image et la surface de I'dbjet dans sa position courante. On 
utilise pour cela des m§thodes inath§matiques de recherche du 
minirnum de la somme des carres de ces hearts par rapport a\ax six 
paraitietres (trois degr6s de libertg en rotation et trois degrSs 

15 de liberty en translaticm) dSfinissant la transformation rigide 
recherchSe. On obtient xane configuration telle qu'illustrge par 
la figure 3B dans laquelle le nodSle 21 « a una position et 
orientation correcte tnsLis une forme encore iraparfaite. C'est 
pourquoi certain rayons r« ne sont pas tangents tnais traversent 

20 le module. 

De preference, la position resultant de la transfonna- 
tion rigide (rotation, translation) est cbtenue par une tnethode 
(algorithme) dite du point le plus proche (iterative Closest. 
Point, ICP) • IM exeinple d'une telle mithode est d^crit dans I'ar- 

25 tide "A Method for Registration of 3-D Shapes" de Paul J. Best 
et Neil D. McKay, paru dans IEEE TRANSACTIONS ON PATTERN ANALYSIS 
AND MACHINE INTELLIGENCE, vol. 14, NO. 2, en f^vxier 1992, et 
dont le contenu est incorporS ici par r6f6rence. 

De la figure 3B k la figure 3C, la transformation est 

30 "non- rigide" (glastique) et consiste a deformer le module 21' 
sans changer son orientation jusqu'a dbtenir 1* image 22 en trois 
dimensions pour laquelle tous les rayons de rStrppro j ectirai r" 
sont a une distance minimale des contours (idScG-^nent tangents) . 
Cette transformation non- rigide est effectu6e, par es^tple et 

35 selon un premier mode de realisation, en calculant les param&tres 
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de deformation a partir d'xon algorithms connu sous la dencxnina- 
tion anglo-saxonne "Down Hill Simplex". Get algorithms a 6t6 
decrit, par exenple, en 1965 par J. A. Nelder et R. Mead dans 
CSOMPXJIER JOURNAL, vol. 7, pages 308-313, dont le contenu est 
5 incorpor§ ici pcur r6f6rence. 

Selon ion mode de realisation pref6r6 de 1' invention, 
I'algorithme utilisS pour determiner la transformation non-rigide 
du modele statistique est bas€ sur la mSthode de Levenberg- 
Marquardt. Le principe de cette methode est de determiner le 

10 minimum d'une fonction non-linSaire multidimensionnelle en utili- 
sant les d^rivges partielles de la fonction par rapport axox para- 
metres de deformation du modele. Cette methode a ete decrite, par 
esffittple, en 1963 par D.W. Marquardt dans JOURNAL OF THE SOCIEIY 
FOR INDUSTRIAL AND APPLIED MATHEMATICS, vol. 11, pages 431-441, 

15 dont le contenu est incorpore ici par reference. La fonction dont 
on cherche le minimum est ici la scinrae des carres des distances 
entre un ensemble de rayons de retrqproj ection s'appuyant sur les 
points de contours et la surface du modele. Les parametres de la 
recherche du minimum sont les coefficients que I'on applique k 

20 chaque mode de deformation du mcxiele statistique, ainsi que cela 
est presente dans 1« article "Building a Ccnplete Surface Model 
from Sparse Data Using Statistical Shape Models: Application to 
Conputer Assisted iQiee Surgery? de Markus. Eleute et Stephane 

Lavaliee, dej^ cite. 

25 Cette methode est appliquee et repetee a chaque 

position approchee jusqu't ce qu'il ne soit plus possible de 
minimiser les ecarts. On enchaine ainsi lan recalage rigide, un 
recalage eiastique, un recalage rigide, im recalage eiastique, 
etc. Une fois la forme en trois dimensions reconstituee, sa 

30 taille est bien entendu egaO-ement connue dans la mesure ou toutes 
les vues en dexxx dimensions ont des tallies connues, de mime que 
la position, dans le referent iel de reference, du modele par rap- 
port aux sources. Le nombre de passes necessaire depend essen- 
tiellement du nombre de parametres caracteristiques utilises pour 

35 mcxieiiser la deformation. 
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On notera que I'ordre dans lequel scait effectuies les 
deux Stapes ci-dessus est inportant. En effet, si can ccraraence par 
effectuer une difortnation non-rigide, on va faire correspondre un 
tnodSle avec une mauvaise orientation qu'il sera tres difficile t 
5 r6cup§rer en raison de la difortnaticxi de forme subie. Toutefois, 
la tninitnisation des hearts peut §tre effectu^e glc±)alenient en 
prenant en cornpte les paramitres rigides et les paramitres de 
dSformations en m&ne temps H chaque iteration du processus mathS- 
matique de minimisation. 

10 On notera que, Sl chaque 6tape iterative du proc6de dans 

laquelle on calcule les hearts entre les contours du meddle et 
les rayons de retroprojection, de nouveaixsc points entrant dans la 
mesure sont d^terminSs cQt€ modele et/ou c6t6 rayons de projec- 
tion. Ceci constitue une distinction importante par rapport a des 

15 proc6d€s connus de reconstitution d' images en trois dimensions 
dans lesquels les points de mesvire sont les m&nes pour toutes les 
iterations . 

une autre caract6ristique de la pr^sente invention est 
que les Evaluations des hearts minimaux ne sont pas nScessaire- 

20 ment faites de fagon ejchaustive pour tous les €ldnents de surface 
du module. En effet, cela conduirait parfois a des calculs trop 
dispendieux m§me si le nonibre d'61§ments de surface peut §tre 
minimis^ en utilisant des Pigments triangulaires; Selon l' inven- 
tion, les gvalioations d 'hearts sont faites par report ^ des 

25 lignes caract&ristiques const ituant des gin^rateurs de contours 
du module. Cela minimise le notribre de points de mesure. Par exem- 
ple, un OS de type f&nur, d^fini par environ 5000 bords de surfa- 
ces triangulaires, est dSfini par environ 300 g6n6rateurs de 
contours . 

30 De preference, 1' analyse effectuEe par le procide de 

1' invention ne touche que des g6nerateurs de contour de 1' image 
en trois dimensions. Ainsi, la surface du module statistique est, 
de preference, fomnee d' elements triangulaires dont certaines 
aretes definissent des generateurs du contour. Un generateur de 

35 contour est defini par les argtes (bords) des triangles qui, pro- 
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jetes sur un plan, definissent le contour (interne ou externa) . 
Le recours aux generateurs de contours permet de reduire consid§- 
rablement (par exenple, d'au moins xon facteur 10) le nombre de 
points a chercher dans le tnodele pour verifier la correspondance 
5 avec 1" image reconstitute. 

La figure 4 illustre la definition d'xjn g&itrateur de 
contour d'une image (raoddle) en trois dimensions dont la sxirface 
est formge d" Pigments de surface triangulaires tels que dScrits, 
par exemple, dans 1' article "Anatomy-based registration of ct- 

10 scan and intraoperative x-ray images for guiding a sutgical 
robot" de A. Gueziec, paru dans IEEE TRANSACTIONS ON MEDICAL 
IMAGING, 17(5), pages 715-728, en octdbre 1998, et dont le 
contenu est incorporg ici par r6f6rence- 

La figure 4 repr^sente, de fagon sitiplifiSe, deux 

15 triangles 30 et 31 dSfinissant vne portion de surface d'un modele 
en trois dimensions et dont 1' arete commune constitue une genera- 
trice 32 du contour. Le calcul a effectuer pour determiner si xane 
arete constitue ou non xme gen6ratrice du contour consiste a cal- 
culer les angles a et p respectifs entre les normales aiox surfa- 

20 ces des triangles 30 et 31 et les rayons 33 et 34 reliant le 
centre de ces surfaces au centre 35 de projection correspondant a 
la position de la source de rayons X, Si I'un des angles est 
infgrieur a 90^ tandis que 1' autre est superieur a 90o, levir 
arete commune 32 est alors ime gengratrice du contour. 

25 La figure 5 illustre le type de mesure effectuSe pour 

determiner l'6cart e entre un raycoi r de retroprojection et un 
point d'un g&ierateur de contour, Cette figure reprgsente, de 
fagon schematique, tine image I de I'dbjet dans un plan P et une 
forme 21 du modele statist ique placee sur le trajet des rayons de 

30 rStroprojection de 1' image I jusqu'^ la source (non representee) . 

Pour chaque rayon de retroprojection r, on cherche les 
points respectifs 41 et 40 du rayon de retroprojection et d'un 
generateur de contour choisi parmi tous les generateurs, pour 
lequel I'ecart e entre ces points est minimal. 
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Un avantage de la pr^sente invention est qu'elle pernet 
une reconstitution beaucoup plus rapide d'line image en trois 
dimensions par rapport aux techniques connues. 

Un autre avantage de la pr6sente invention est qu'elle 
5 permet vn alignement correct du modele meme dans des zones oQ on 
ne dispose pas de donn6es de projection par les images en deux 
dimensions (par exetrple, certaines courbures internes) . 

On not era que le contour de surface pourra, le cas 
6ch6ant, etre affine au moyen d'un palpeur mecanique, qptique ou 
10 magngtique foumissant des coordonn^es spatiales dans le m§me 
r^ferentiel que les images en deux dimensions. On obtient alors 
des points supplementaires qui peuvent §tre utilis§s dans la 
recherche de 1' image en trois dimensions. Le recours §. un palpeur 
peut servir, par exenple, a diminuer le nowibre de vues en deux 
15 dimensions nScessaires en donnant une information en trois dimen- 
sions. 

On notera Sgaloment que le proc6d6 de 1» invention peut 
s'appliquer §i plusieurs contours surfaciques en trois dimensions 
constituant im ou plusieurs cbjets, Par exancple, on peut ainsi 

20 reconstituer simultanSnent la surface exteme et la surface cor- 
ticale interne d*\m os. On peut 6galCTient reconstituer la forme 
de plusiexirs os participant Si un articulation et visibles simul- 
tan6nent sur les images radiologiques (par ex^ple , le tibia et 
le f &nur si on s ■ intSresse au genou) , ou encore la forme de plu- 

25 sieurs fragments osseux d'un meme os. Le moddle recherche doit, 
dans ce dernier cas, contenir les transformations rigides entre 
chaque contour surfacique en trois dimensions. Dans d' autre cas 
d' application a plusieurs cbjets, il faut bien entendu connaitre 
les dif ferents modeles statistiques et les ^entuelle transforma- 

30 tions rigides ou Slastiques entre eux. 

Bien entendu, la pr6sente invention est susceptible de 
diverses variantes et modifications qui apparaitront iL I'hcmme de 
I'art. En particulier, le nombre de vues en deux dimensions 
devant §tre utilisS pour la reconstitution depend de la precision 

35 souhait§e et de la cooplexitS du moddle. Dans certeiins cas, xine 
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seule vue en deux dimensions pourra suffire. De plus, bien que 
1' invention ait §te decrite en relation avec une source de rayons 
X tnobile, celle-ci pourra etre renplacee par plusieurs sources 
fixes, pourvu de respecter la contrainte d'dbtenir des images en 
5 deux dimensions dans un meme rSferentiel. En outre, la mise en 
oeuvre de 1' invention fait bien entendu appel a des techniques de 
traitement niamerique d* images dont la recilisatian pratique est ^ 
la port^e de I'homtne du metier a partir des indications fonction- 
nelles donnees ci-dessus et dans las publications cities en refe- 
10 rences. 




19 

REVENDICATIONS 

1. Proc6cl6 de reconstitution d'une image en trois 
dimensions repr^sentant les contours surfaciques d'au tnoins un 
objet (1), k partir d'au moins une vue en detix dimensions de cet 
dbjet prise aux rayons X, caracterise en ce qu'il consiste a : 

5 determiner la position de la source (7) de prise de 

vues dans vn riferentiel de r^f^rence ; 

s§lectionner un moddle predefini constituant une forme 
moyenne (21) de 1' dbjet / et 

de fa<?on iterative jusqu'^ ce que les contours du 
10 modele soient tels que les 6carts, entre des rayons de rStropro- 
jection des contours de 1' image en deux dimensions depuis la 
source et la sxirface du moddle, soient minimaux : 

sSlectionner une orientation et une positicxi du module 
dans le r^fSrentiel de r^fSrence, puis 
15 sSlectionner une dSforroation du raodile pour modifier 

ses contours en trois dimensions. 

2. Proc6d6 selon la revendication 1, caract6ris§ en ce 
que le mod&le est dbtenu a partir d'une population d'dbjets pour 
laquelle on recherche la correspcxidance statistique ccmraune S. 

20 tous les objets pour determiner une forme moyenne et les deforma- 
tions principales par rapport a cette forme moyeaine, de fagon 4 
disposer d'au moins un moddle statistique. 

3. Proc6d§ selon la revendication 2, caracterise en ce 
que les Stapes de selection iteratives consistent a faire subir 

25 au moddle statistique, successivement, une transformation rigide 
modifiant sa position et/ou son orientation et une deformation 
non-rigide modifiant ses contours surfaciques. 

4. Precede selon I'une quelconque des revendications 1 
a 3, caracterise en ce que les contours de 1' image en deux dimen- 

30 sions sont obtenus automatiquement en projetant le modSle dans le 
plan de 1' image en deux dimensions, et en deformant les contours 
projetes de fa^on a les faire coincider avec les points de fort 
gradient en niveaux de gris de 1' image en deux dimensions- 
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5. Precede selon la revendication 4, caracterise en ce 
que la determination autcxtiatique des contours de I'irnage en deux 
dimensions est effectuee de fa?on iterative, chaque iteration 
gtant intercalee entre deux iterations successives des Stapes de 

5 selection. 

6. Proc6d§ selon I'une quelconque des revendications 1 
a 5, caracterise en ce qu'il consiste S. determiner, dans le refe- 
rentiel de reference, des coordonnSes en trois dimensions de 
points de I'objet, de fa^on k disposer de points de reference 

10 suppiementaires pour les etapes de selection it6ratives de posi- 
tion, orientation, et deformation. 

7. ProcedS selon I'xane quelconque des revendications 1 
a 6, caracterise en ce qu'il consiste a utiliser plusieurs images 
en deux dimensions pour lesquelles les positions respectives de 

15 la source de prise de vues sont toutes determinees dans le refe- 
rent iel de reference, et i effectuer les etapes de selections 
iteratives en tenant cornpte des rayons de retroprojection des 
contours de toutes les images en deux dimensions. 

8. Procede selon la revendication 7, caracterise en ce 
20 que le norabre d' images utilise est fonction de la precision 

souhaitee . 

9. Procede selon I'une quelconque des revendications 1 
a 8, caracterise en ce que la surface du module est constituee 
d' elements de triangle, lesdits ecarts etant mesxires par rapport 

25 §i des points de certaines argtes const ituant des generateurs du 
ccffitour en trois dimensions. 

10. Precede selon I'une quelconque des revendications 1 
a 9, caracterise en ce qu'il est applique k la reconstitution des 
contours surfaciques de plusieurs dbjets lies entre exjx par des 

30 relations de transformation rigide et/ou eiastique. 

11. Procede selon I'une quelconque des revendications 1 
a 10, caracterise en ce qu'il est applique k la reconstitution 
d'images d'os. 




12. Systdme de traitement d'itnages, caractSris^ en ce 
qu'il cottporte des inoyens pour la mise en oeuvxe du precede selon 
I'une quelconque des revendications 1 ^ 11. 
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Fig 3A 




7 Fig 3B 
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Fig 5 



